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ABSTRACT 

Using a Newtonian approximation we developed a quantitative criterion for the col- 
lapse of a spherical distribution of matter under an isolated texture field. In particular, 
we found the evolution of an overdense region is strongly determined by two param- 
eters: the energy scale of symmetry breaking (ry) and the initial radius of the system. 
Applying our collapse criterion to typical galaxy scales we verified the formation of 
10"Mq objects at z < 9 and 10^'^Mq objects at z < 5. 

Key words: Cosmology: structure formation; cosmic fields: global textures. 



1 INTRODUCTION 

The emergence of structures in the universe is one of the 
most important problems of modern cosmology. The diffi- 
culty is to understand how an initially homogeneous mass 
distribution evolved into its present clumpy state. Gener- 
ally, it is accepted that structures were initiated by small 
density fluctuations formed at the early universe, and that 
the subsequent clustering was produced by gravitational in- 
stability (e.g. Kolb & Turner 1990). At the moment, there 
are two theoretical approaches proposed to explain the ori- 
gin of the primordial inhomogeneities. In the context of in- 
flationary models, they result of quantum-induced Gaussian 
fluctuations amplified over an accelerated phase of cosmic 
expansion. On the contrary, in topological defect scenarios, 
non-Gaussian fluctuations are produced during spontaneous 
symmetry breaking at phase transitions in the early universe 
(e.g. Vilenkin & Shellard 1994). 

Several types of topological defects can be formed de- 
pending on the kind of symmetry that is broken. Here, we are 
interested in a speciflc topological defect formed during the 
symmetry breaking of non-Abelian groups: a global texture. 
This defect has been proposed as a possible source of galaxy 
formation in the universe (Turok 1989). The development 
of structures based on this topological remnant is related 
to the evolution of textures knots. On scales larger than 
horizon, textures knots are regions where the scalar field 
winds around the vacuum manifold in a non-trivial way. As 
knots come inside the horizon, they collapse at the speed of 
light down to an infinitesimal scale where they unwind them- 
selves emitting spherical waves of outgoing radiation. This 
collapse produces an overdense region onto which matter is 
attracted, possibly leading to the formation of non-linear 
objects (Turok & Spergel 1990). 



The most desirable way to study structure formation 
produced by textures knots is through numerical simula- 
tions like those carried out by Park, Spergel & Turok (1989), 
Spergel et al. (1991) and Gen et al. (1991). AU these stud- 
ies indicate important effects of textures on the distribution 
of matter, like coherent large-scale flows and clustering of 
galaxies fairly consistent with the correlation function ob- 
served in the nearby universe. However, we should keep in 
mind that the dynamics of textures is highly non-linear since 
the decoupling era (Albrecht, Battye & Robinson 1997), and 
that it is not simple to trace the evolution of the fluctuations 
over this long time interval. Thus, all these numerical works 
are strongly dependent on the choice of the initial conditions 
for the N-body simulations (Deruelle et al. 1997). 

Actually, this is not the only problem with the tex- 
ture scenario. The work of Pen, Seljak & Turok (1997) 
shows that the texture model has strongly suppressed acous- 
tic peaks in comparison to current observations of GMB 
anisotropics. Also, when COBE normalized, the model pro- 
duces a low value of the rms density variation on scales of 
8 Mpc (erg ~ 0.23), while a higher value of this parameter 
(closer to unity) is necessary in order to get agreement with 
the observed galaxy clustering. These results clearly pose se- 
rious difficulties to the texture scenario. However, this does 
not mean that the model should be completely discarded at 
the present time, since a decisive, model discriminating test 
will be possible only after the analysis of the high resolution 
GMB maps which will be available in the next years from 
the MAP and PLANGK satellite missions. Hence, for the 
moment, it remains valid to work with texture knots as a 
viable way to form galaxies in the universe. 

Indeed, an alternative way to study structure formation 
via global textures is to develop simple analytic models to 
describe the general properties of the defect dynamics. For 
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instance, Gooding, Spergel & Turok (1991) studied the for- 
mation of bound objects induced by a number of collapsing 
knots in a Qo = 1 CDM-dominated universe. An important 
result of this study is the prediction of early spheroidal for- 
mation in the universe, suggesting that by 2 ~ 50, ~3% of 
the mass of the universe has formed non-linear objects with 
mass greater than 10® Mq and that most objects larger than 
10^^ Mq formed by z ~ 2 - 3. 

In the present work, we use a Newtonian approximation 
to investigate the formation of self-bound gravitational sys- 
tems seeded by global textures. However, instead of applying 
a statistical approach to the contribution of many individ- 
ual knots, we are only interested in the local influence of a 
single texture knot on the process of galaxy formation. Our 
aim is just to obtain a quantitative criterion for the collapse 
of a spherical distribution of matter under a texture field. 



2 THE COLLAPSE CRITERION 

Distributed field gradients induced by a number of global 
textures are supposed to give a larger contribution to den- 
sity inhomogeneities than that produced by single isolated 
knots (e.g. Vilenkin & Shellard 1994). If this is correct, then, 
the effort to probe the process of structure formation seeded 
by these defects will require large dynamic range simulations 
in computers of latest technology. However, the demand for 
such numerical works does not mean that analytical devel- 
opments are unnecessary. On the contrary, they are particu- 
larly useful to give a supplementary and intuitive picture of 
the defect dynamical contribution on structure formation. 
Specifically, in this work, we study the process of accretion 
of matter onto a single texture knot. This problem has been 
studied by several authors over the years. But now, our em- 
phasis is on the Newtonian treatment and the possibility of 
taking a region where the assumption of an isolated texture 
is reasonable. 

In order to do that, let's consider an initially spheri- 
cal and homogeneous distribution of matter with density pt 
when the texture field is turned on. The Newtonian density 
associated to the texture can be obtained by using the Ein- 
stein equations in the weak field approximation. By numer- 
ically solving the Barriola-Vaschaspati equations for a self- 
gravitating texture, Gueron & Letelier (1997) found that 
the weak field condition is fulfilled and the corresponding 
Newtonian density is 

where 77 is the energy scale of symmetry breaking. This den- 
sity is negligible except when the texture size, t, is small 
compared to the physical scale r. At the same time, the re- 
gion where is reasonable to take account a single texture 
configuration, ignoring the effects from other textures, is 
limitted by 



{At^ + r^^^ <, H-^ 



(2) 



(Notzold 1991), where At — t, — t, t» is the time at which 
the texture is collapsed and H is the cosmic expansion rate. 
For simplicity, we will only study the behaviour of regions 
where r ;» |i|, so we can use an effective density defined as 



For such regions, |Af| < r and so condition (2) is simply 
r ^ H~^ l\pl. The total density within r is p(r) = pt + Pb 
and the density contrast can be expressed by 



(4) 



5p ^ p{r) - Pb ^ PT 
P Pb Ph 

so that the motion of a thin shell of particles located at r is 
governed by the equation 



d> 



GM 



, where M 



iTvr' 



dt^ r2 ' 3 

and the average density contrast is 



Pb{l + S) 



4y^j.3 



4:Tvr dr. 



(5) 



(6) 



As usual, it will be assumed here that the mass within 
the shell is constant. In fact, this is strictly correct only after 
a time interval of ~ i, , when the texture mass is converted in 
massless goldstone bosons, leaving behind the other material 
components: baryonic and dark matter. Thus, we are tacitly 
assuming that t, is much smaller than the dynamical time 
of the system, but large enough to allow the texture mass to 
dominate the first steps of the gravitational instability. In 
idealizing the process this way, we hope to make it well-posed 
mathematically, although a justification of the idealization 
will be only possible from the agreement (or not) of our 
predictions with the data. 

Now, supposing that the peculiar velocities at t = ii 
are negligible, the first integral of Eq.(5) will give the total 
energy of the shell and can be written as 



E = KiQilQ^^ - {l + 5i) 



(7) 



(e.g. Padmanabhan 1993). The subscript i indicates the ini- 
tial time ti when the radius of the shell is ri, the kinetic 
energy is K^ and the density parameter is f2i. Naturally, the 
condition E < for a self-bound system is 



Si > {nr^ -l) = 3c. 

For a dust configuration (p = 0), we have 

!^o(l + Zi) 



Qi 



1 + ZiOo 



and so 



5c = 



fio 



no{l + Zi) 



(8) 



(9) 



(10) 



(e.g. Borner 1988). Therefore, for the case of a flat universe, 
5c — and any overdense region with Si > will collapse. 
On the other hand, for the case of an open universe, only 
regions where Si > Sc can produce self-bound gravitational 
systems. In this case, introducing definition (3) in Eq. (6) we 
directly find 



Si = 



24r?2 



Pb{ti)rf 



(11) 



Now, it is easy to foresee the destiny of a spherical dis- 
tribution of matter under a texture field when fio < 1- In- 
troducing Eq.(ll) in condition (8), it is simple to show that 



n < 



24772 t^i 



Pt = 8rj /r 



(3) 



Pb(ii)(i-no' 

Recalling that 



(12) 
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Pb(ii) = 0,ipc{ti) 
we reach 

Sri I Gtt 
ri < 



HiV 1 - r^i 



and 



Re 
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(14) 



where R^ is the critical radius for the collapse. Once more, 
assuming a Friedmann model with vanishing pressure and 
with no cosmological constant 

Hf = J/o [^0(1 + Zif + (1 - no)(l + Zif] (15) 

(e.g. Sandage 1961), we can finally express the collapse cri- 
terion for any overdense region under a texture field at any 
epoch zi: 



n < 



GTT{l + Zino) 



8^ 

Ho V (1 - f^o)[no(l + zi)3 + (1 - no)(l + zi)2] ■ 



(16) 



3 THE CHOICE OF THE ENERGY SCALE 

From Eq.(16) we see that, assuming specific values for 
Ho and f^o, the critical radius for the collapse will depend 
basically on the the energy scale of symmetry breaking. Gen- 
erally, global textures arise after a phase transition during 
the break of symmetry described by the Grand Unified The- 
ory (GUT), when the strong interaction is unified with the 
eletroweak force. Current estimates of this unification lead 
to an energy scale of ~ lO^'^ GeV (Gleiser 1998). In the same 
way, if the observed fluctuations in the cosmic background 
radiation are due to textures, we must have {5T/T) ~ SGr/^ 
which implies rj ~ 10^'' GeV in order to obtain an amplitude 



5T/T 



10" 



consistent with the COBE measurements, 
(e.g. Vilenkin & Shellard 1994). 

In view of these arguments, the energy scale of 10^® 
GeV seems to be the natural choice of 77. However, let's 
proceed here to find a new and independent argument for 
this choice. First of all, assuming fio ~ 0.2 and h = 0.75 
(h = _ffo/100 km s~^ Mpc~^), we plot the behaviour of _Rc 
as a function of the epoch Zi using r; as a free parameter 
(see Figure 1). Note that for earlier epochs Re is smaller 
(independent of 77) indicating that the first bound objects 
to be formed via textures knots would have small sizes and 
masses. In fact, isocurvature fiuctuations like those produced 
by topological defects can evolve only after the decoupling 
(z ~ 10^), when the Jeans mass is ~ W^Mq. Thus, clumps 
with this mass should form first and evolve into larger sys- 
tems through gravity. 

Obviously, the details of the evolution of clumps up to 
form typical galaxies and clusters is a rather complex pro- 
cess, involving the thermal history of the gas, star formation 
phenomena and the mutual interaction between the clumps. 
It is beyond our aim to address this problem as a whole. 
Instead, we proceed to simplify it by supposing that a pro- 
togalaxy just consist of a regular distribution of matter (pos- 
sibly sub-galactic clumps) under a texture field where we can 
apply our collapse criterion. 

That being the case, an important quantity to be found 
is the mass Mc related to a given radius Rc{zi). For a home- 
geneous universe, the non-relativistic mass within a physical 
scale A is given by 



Figure 1. Plot of _Rc(Mpc) as a function of the redshift zi and 
the energy scale of symmetry breaking rj, where we have used 
h = 0.75 and Co = 0.2. The dotted line corresponds to the limit. 



M(A) = 1.45 X lo"(J7o/i )AmpcM0 



(17) 



Figure 2. Plot of M (in solar units) as a function of the redshift 
Zi and the energy scale of symmetry breaking rj. The dashed line 
indicates the behaviour of the Jeans mass with z; . We have used 
h = 0.75 and Ho = 0.2. 



(e.g. Padmanabhan 1993). This quantity is conserved during 
the expansion and corresponds to the uncoUapsed mass of a 
physical scale A. If we take A = _Rc, it is possible to study 
the behaviour of Mc as a function of the epoch Zi and the 
parameter rj. In Figure 2 we see that, for a same epoch Zi, 
the value of A/ will strongly depend on 77. From the com- 
parison between Mc and the Jeans mass, we also note that 
the lowest energy scale which is able to form self-bound sys- 
tems over all the astrophysical range (10^ — IO^^AIq) at any 
pos-recombination time Zi is rj = 10^^ GeV. At this energy 
scale. Re is well inside H'^^ 1^/2 (see Figure 1), which defines 
regions where we can consider the effects of isolated texture 
knots. This result can be taken as a qualitative argument 
for the choice of 77 = lO'^^ GeV and so we are assuming this 
value for the rest of the work. 
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Figure 3. Plot Si as a function of Zi. Tlie dotted line is the limit 
2.78(1 + 2i)'^ 



4 GALAXY FORMATION 

In order to apply our collapse criterion to the specific 
problem of galaxy formation, we should follow the evolu- 
tion of overdense regions containing typical galaxy masses 
(10" -10^^ Mq). If these respect the condition (8), they wiU 
expand to a maximum radius, then collapse and eventually 
virialize. The dynamical properties of the resultant systems 
can be estimated from the spherical model for non-linear 
collapse. Thus, at the epoch 



(l + 2col) =0.365i(l + Zi) 



(18) 



the dissipative component of the system (i.e. baryons) 
reaches the virialization with the following approximated 
properties: 

-2/3, 



163(1 + zc 






and 



126(1 



^l/2^2/3f 



(M/IO^^Mq)^/^/!" 



1/3,1/3 



z,oi)-''nt,'-'{M/WMQy"h 



^kpc (19) 



kms"^ (20) 



(e.g. Padmanabhan 1993). 

Before applying these expressions to galaxy scales, it 
would be interesting to assure that the initial density con- 
trast Si can form a self-bound system before the present 
epoch, Zcoi > 0. We directly find such a condition from 
Eq.(18): 

Si > 2.78(1 + Ziy'^ . (21) 

For the case of a flat model this is a simple condition for 
the collapse. However, for an open universe, we also should 
have Si > 5c in order to collapse the overdense region. Using 
Eq.(lO) and taking the limitting case when Sc — 2.78(1 + 
Zi)~^ , we conclude that f2o < 0.26 is the necessary condition 
for the collapse with Zcoi > 0. 

In Figure 3, we plot Si as a function of Zi for initial radii 
that correspond to typical galaxy masses. The dotted line 
is the limit 2.78(1 -I- Zi)~^ . Note that for M = lO^M©, the 
collapse is only possible for Zi ^ 33 (zcoi ^ 9) if Slo = 0.2 
and Zi <, 25 (^^oi ^ 8) if ilo = 1-0; while for M = lO^^Af©, 
Zi <> 15 {zcoi ^ 5) if Qo = 0.2 and Zi <, 11 (zcoi ^ 4)if 



Figure 4. cr — r plane of the virialized systems. The boxes delimit 
typical values found in the literature for spiral (hatched box) and 
elliptical galaxies. Curves A and B refer to 10^^ Mq objects and 
curves C and D refer to 10^^ Mq objects. 



Slo = 1.0. These ranges in Zi produce associated loci in the 
a — r plane (see Figure 4) . These loci show that our simple 
model texture knot-|-spherical collapse is able to produce 
self-bound gravitational systems with dynamical properties 
quite similar to real galaxies. In fact, the comparison with 
data compiled from the literature shows a significant agree- 
ment between our predictions and typical galaxies, while it 
is also clear in Figure 4 an important trend towards larger 
radii and lower velocity dispersions in our bound systems. At 
the same time, the epoch when such objects are presumed 
to reach the virial equilibrium (zcoi) is consistent with the 
results of several studies which point out an epoch of galeixy 
formation hy z ^ 5 (e.g. Peebles 1993; Pahre 1998). 

Now, by extending the approach developed in this sec- 
tion to a wider astrophysical range (10^ — 10^'' Mq), it is 
possible to depict an intuitive picture of the structure forma- 
tion process when it is seeded by textures knots. The results 
for both open and flat universes are summarized in Tables 
1 and 2, respectively, where the columns are as follows: (1) 
is the mass of the systems in solar units; (2) is the initial 
redshift we take; (3) is the initial radius of the region to be 
collapsed in Mpc; (4) is the initial average density contrast; 

(5) is the redshif at which the systems reach the virialization; 

(6) is the time interval of the collapse and virialization pro- 
cess (the dynamical time) ; flnally, (7) and (8) are the radius 
in kpc and the velocity dispersion in km/s of the virialized 
systems. 

A simple analysis of the data assembled in these tables 
reveals some points we should emphasize. First of all, note 
that structure formation ceases at earlier times in the low 
Qq universe in comparison to the flat case. This is expected 
as long Eis matter has a smooth distribution in high f2o uni- 
verses so that it contributes to the expansiong rate but does 
not contribute to the gravitational potential and so it slows 
down the growth of fluctuations. At the same time, note that 
the process of structure formation in our model begins at a 
much earlier epoch in comparison to standard CDM models, 
where most of the objects with 10^ A^q are formed at the 
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epoch z = 50/fe (6 is the biasing factor) (e.g. Peebles 1993). 
Therefore, even an unbiased CDM model would present a 
less advanced structure formation stage for a given redshift 
in comparison to our model. 

Such an early self-bound system formation may have 
serious consequences. In particular, we form the first gener- 
ation of 10*^ Mrs, objects at z = 135 {flo = 0.2) and z = 112 
(Qq = 1.0). These objects could be related to an early star 
formation phase in the universe and its consequent reioniza- 
tion. At the same time, these early spheroids may develop 
massive black holes onto which matter would be accreted, 
possibly forming an early generation of active galaxies. 

Still concerning to the chronology of structure forma- 
tion, Gooding, Spergel & Turok (1991) found that by z ~ 50, 
~ 3% of the mass of the universe has formed non-linear ob- 
jects of mass greater than lO'' Mq and that most of ob- 
jects larger than 10^^ Mq form by « ~ 2 — 3. These results 
should be compared to our findings, while in the present 
work we have no statistics on the fraction of matter which 
is collected in bound objects of a specific mass in a given 
redshift. Note, however, that by z ~ 50 we already formed 
systems of lO'' — 10* Mq and by z ^ 3 we form objects 
with 10^^ — 10^'* Mq. Thus, in some measure, our results 
are consistent with Gooding, Spergel & Turok (1991) 

On the other hand, a disadvantage of our model refers 
to the velocity dispersions of the collapsed systems which 
are sistematically lower than those observed in astrophysi- 
cal objects of same mass. Up to galaxy scales this effect is 
accompanied with larger radii and so it can be taken as a 
simple consequence of the virial theorem. On the contrary, 
the lower velocity dispersion effect seems to be particularly 
serious in larger scales, like galaxy clusters, where the model 
predicts objects with radii similar to those found in real 
clusters. Actually, this is the opposite one could expect in 
a texture scenario, where the positive skewness of the mass 
fluctuations induces higher velocity dispersions for objects 
of a given mass than (as an example) in standard CDM 
models (e.g. Bartlett, Gooding & Spergel 1993). We think 
this trend is probably due to the fact we are considering 
only the dissipative baryonic component of matter inside 
rvir. By adding dark matter within this radius we could 
recover values of velocity dispersion nearer to reality. For 
instance, in the case of a galaxy cluster scale (10^* A'^©), if 
the total mass (dissipative plus non-dissipative components) 
inside rvir reached 10^^ Mq the velocity dispersion of the col- 
lapsed system would be 583 km/s (f2o = 0.2) or 1475 km/s 
(fio ~ 1.0), which could be taken as more "normal" val- 
ues. However, if this is the right solution to the problem, we 
would have to assume specific distributions for the luminous 
and dark matter components when considering different ob- 
jects. At the same time, we cannot discard the possibility 
that the simplicity of our model texture knot-|-spherical col- 
lapse is not taking account all of the physics relevant to the 
structure formation process. Indeed, we are ignoring interac- 
tions between clumps, star formation and hydrodynamical 
effects. These physical mechanisms are probably important 
over the systems evolution and their inclusion would change 
(to some extent) the outcomes of the model. 



Table 1. Properties of self-bound gravitational systems seeded 
by textures knots in a Qq = 0.2 universe. 



M 
10'^ 



<5i 



-^col 



tr 



1000 0.04 0.0067 135 7.7 x 10" 



1.4 



4.6 



10^ 


700 


0.09 


0.0038 


80 


1.6 X 10^ 


7.3 


7.5 


10« 


350 


0.18 


0.0076 


50 


3.3 X lO'^ 


16 


13 


10» 


150 


0.39 


0.0205 


27 


7.6 X lO'^ 


34 


21 


lOio 


70 


0.85 


0.0416 


16 


1.5 X 10* 


73 


35 


lOii 


32 


1.84 


0.0886 


9 


3.0 X 10» 


155 


58 


1012 


14 


3.87 


0.2131 


5 


5.4 X 10« 


333 


98 


10" 


6 


8.57 


0.4273 


2.7 


8.3 X 10® 


680 


163 


10" 


2.4 


18.4 


0.8100 


1.2 


1.2 X lO'^ 


1465 


271 



Table 2. Properties of self-bound gravitational systems seeded 
by textures knots in a Qq = 1-0 universe. 



M 
10^ 



<5i 



■^col 



''col 



1000 0.02 0.0054 112 6.9 x 10® 



1.4 



12 



|7 


550 


0.05 


0.0052 


66 


1.5 X lO'' 


4.3 


20 


|8 


260 


0.11 


0.0100 


39 


3.2 X 10^ 


9.2 


34 


|9 


115 


0.23 


0.0261 


23 


6.7 X 10^ 


20 


56 


|10 


55 


0.50 


0.0491 


13 


1.4 X 10« 


42 


93 


|11 


25 


1.07 


0.1070 


8 


2.5 X 10** 


91 


157 


|12 


11 


2.31 


0.2330 


4 


5.7 X 10« 


190 


261 


|13 


5 


4.98 


0.4020 


2 


1.1 X lO'' 


425 


422 


|14 


1.6 


10.7 


1.0000 


0.7 


8.5 X 10* 


857 


685 



5 DISCUSSION AND SUMMARY 

In this work, we developed a quantitative criterion for the 
collapse of a spherical distribution of matter under a sin- 
gle texture field. Further, we applied this criterion to form 
galaxies in an idealized smooth universe where texture knots 
can be taken as isolated sources for gravitational collapse. 
Despite the simplicity of the model, the properties of the 
systems which collapse and virialize respecting our criterion 
are similar to real astrophysical objects. A drawback of the 
model is the trend towards larger radii and lower velocity 
dispersions in our bound systems in comparison to typical 
values found in the literature for real objects of same mass. 
This trend is specifically significant for velocity dispersions 
and could be related to the quantity of dark matter inside 
the virialization radius. However, we would like to empha- 
size that the most important achievement in this work is 
the determination of a collapse criterion using simple argu- 
ments and the possibility of applying it in a direct way to 
predict astrophysical objects. This means that local effects 
of isolated knots may be important to the development of 
structures in the universe and suggests that individual tex- 
tures (not only their dynamical collective effect) should be 
taken account in galaxy formation models. 

In the following, we summarize the main results of this 
work: 

(1) We showed that, under specific conditions, the des- 
tiny of an overdense region under a single texture field will 
basically depend on its initial size and the energy scale of 
symmetry breaking rj (if f2o < 1) . For those regions smaller 
than 7?c, the system will expand to a maximum radius and 
then inevitably collapse. 

(2) We found an additional argument for the choice of rj 
as lO'^® GeV. This is simply the lowest energy scale which 
allows the formation of structures on scales relevant to as- 
trophysical objects (10*^ - IO^^Mq) at any time Zi < 1000. 
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(3) We showed the loci of objects with masses 10^^ — 
10^^ Mq in the a — r plane agree well with the observed 
properties of real galaxies. 

(4) Our results also indicate an epoch of galaxy formation 
fairly consistent with other independent studies which point 
out it should have been by z ~ 5. 

(5) Our model produces a chronology of structure forma- 
tion which begins at earlier times in comparison to other 
models, in particular CDM models. 

(6) The existence of spheroids by 2: ~ 100 suggests the 
possibility of associated phenomena like early star forma- 
tion process, reionization of the universe, production of black 
holes and active galaxies. 
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